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1. Introduction 

In e + e~-annihilation the total hadronic cross section is given by the 
probability to produce an initial qq pair, which is determined by QED. 
The effects of the strong interaction is here only a correction with relative 
magnitude a s jix. The emission of a gluon cascade is, however, essential 
for the properties of the final state, although it does not change the total 
cross section. Here angular ordering is crucial for the result, and the dipole 
formulation particularly convenient. 

DIS and hadronic collisions are more complicated. In a high energy ep 
or pp collision the initial partons in a target proton develop virtual parton 
cascades. A projectile can interact with any of the partons in the cascade, 
which implies that the total cross section grows with increasing collision 
energy. The development of the cascade in this "initial state radiation" 
determines the inclusive total and elastic cross sections, but for exclusive 
final states also "final state radiation" has to be added. In the initial state 
radiation the virtualities are spacelike. The final state radiation is more 
similar to the cascades in e + e~-ann.; it does not change the inclusive cross 
sections, and the virtualities are timelike. Thus in DIS and hadronic colli- 
sions we have two different problems, the total cross section and the final 
state properties. There are also two different hard scales, Q 2 and s, while 
in e + e~-ann. there is only one, Q 2 = s. 
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At high energies and small x, gluon cascades and the 1/z pole in the 
splitting function are most important. For large Q 2 this leads to the DLL 
approximation, and for limited Q 2 to BFKL evolution. 

The high density of partons in a proton also implies that at high energies 
the projectile may interact with more than one parton in the target. As the 
total interaction probability must not be larger than one, the effective gluon 
density must "saturate" . At high energy the impact parameter is related to 
the conserved angular momentum, b ~ L/k. The interaction probability for 
fixed b is therefore limited by 1, and a description of multiple interactions 
and saturation is most easy in impact parameter space. 

The outline of these notes is first particle production in e + e - -ann. with 
timelike cascades and hadronization in Sees. [2] and [3l followed by small x 
evolution in SecHJand saturation in Sec.[5l and finally a discussion of dipole 
models for high energy pp collisions and DIS in Sec. [6j 



2. Timelike cascades 

2.1. Bremsstrahlung 
In classical electrodynamics the radiation of bremsstrahlung photons is 



given by the expression (see e.g. Ref. pQ 

d 3 k 



d 4 xj(x)A*(x) 



(1) 



For a charged particle moving along the trajectory x = r(i) we get the 
current (the charge is denoted g, as the result is the same in QCD) 

j=gv(t)6(x.-r(t)). (2) 

With a photon field A ~ g e — *(<«>£— kx) we g n( j ) after division and multiplica- 
tion by (1 — nv) and a partial integration, the amplitude 



where 



M = J d 4 xj(x)A*(x) =ig J dt — e 

ev(t) 



iuj(t— nr(t)) 



X 



n = k/cj. 



(3) 



(4) 



uj(1 - nv(t)) 

For soft emissions (small ui) the exponential is approximately constant in 
regions where dv/dt ^ 0, which implies that 

M oc (X f -Xi), (5) 



dn ~ a — dfl 

00 



ev f 



ev, ; 



1 — nvi 



1 - nv,- 



(6) 
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where v, and Vf are the velocities before and after the radiation. (For large 
oj the emission is, however, sensitive to details in the trajectory.) 



2.2. Dipole radiation 

For emission from pair production of a positive and a negative particle, 
moving along trajectories r + and r_, we get the current 

t < : j = 

t>0: j = +0v+(t)5(x-r + (t))-0V_(i)5(x-r_(i)). [t) 

Thus we see that we get the same result as in Eq. (jSJ), only with the re- 
placements Vf — > v + and v, — > v_. 

In the cms system the matrix element for photon emission becomes 



9 • 2 

oj z sin 



9 (p+k)(p-k) 



We note that the last expression is relativistically invariant, and thus can 
be used in any Lorentz frame. 

We can compare this result with the expressions from the relevant Feyn- 
man diagrams. The two factors in the denominator in Eq. ([SD correspond 
to the propagators l/(p+ + k) 2 = l/[2(p + fc)] and l/(p- + k) 2 = l/[2(p_fc)] 
obtained when the photon is emitted from the positive and negative parent 
respectively. Coherent emission from the two parents give the "dipole for- 
mula" in Eq. ([8]). Denoting the particles with momenta p+, p_, and k by 
the numbers 1,2, and 3, and defining Sy = (pi +pj) 2 , we also get (including 
a proper factor 1/V) 



dn= ads 13 ds 23 = ad^ dy 



where 



k 2 ^ y= I ln £23 (1Q) 

S 2 S13 



represent the transverse momentum and the rapidity in the dipole rest 
frame. For gluon emission in QCD we get the same expression with the 
finestructure constant a replaced by N c a s /2. (For radiation from quarks 
there is a suppression factor (1 — 1/N 2 ), which is not present for dipoles 
formed by gluon charges.) 
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2.3. Angular ordering 
With the help of Eq. ([8]) the dipole emission can also be written 

dn ~ a^-d^l Ql2 with = 1 — n^n, = 1 — cos 9a, (11) 
uj a 13 a 23 

where rij is the direction of particle i, and 6ij is the angle between rij and 
xij. As in Eq. ([9]) particle 3 corresponds to the emitted photon or gluon. 
The last factor can be rewritten in the form 

= \[X 1 +X 2 ] (12) 

The first term in the parenthesis (X\) is non-singular when 023 — > 0. Av- 
eraging this term over the azimuth angle, around ni, keeping the polar 
angle #13 fixed, we get 

_L f Xl # = —e(6 12 - 13 ) (13) 

2ir J ai3 

A similar expression is obtained when averaging X 2 for fixed angle #23 • Thus 
approximating X\ and X 2 by these averages, the emission corresponds to 
independent emission from two emitters within the angular ordered regions 
#13 < #12 and #23 < #12 respectively. 

2.4- More gluons 

The emission of two gluons is considerably more complicated. In a com- 
pressed form the lowest order result for qqgg final states takes three full 
pages in Ref. [2j. However, when the emissions are strongly ordered, i.e. 
when P4 <C p 3 -C W, where p 3 and p$ are the gluon momenta, the result 
factorizes, and thus simplifies considerably. In a semiclassical picture the 
hardest gluon is emitted first from the qq dipole. This gluon carries away 
colour charge and thus changes the current responsible for subsequent softer 
emissions. If the first emission produces e.g. a red quark, a blue-antired 
gluon, and an antiblue antiquark, then the red-antired charges radiate co- 
herently as a colour dipole formed by the quark and the gluon. In the 
rest frame of this dipole the distribution is also given by the expression in 
Eq. Q. In the same way the blue and antiblue charges radiate coherently 
as a colour dipole formed by the gluon and the antiquark [3]. (There is 
also a colour-suppressed term corresponding to a dipole spanned between 
the quark and the antiquark, with relative weight —1/N^.) The emission 
of a gluon with transverse momentum k± is determined by an average of 
the current in Eqs. ([H [3]) over the "Landau-Pomeranchuk formation time" 



Q12 1 Q12 ~ Q13 + «23 
«13G23 2 _ 013023 



+ (1 o 2) 
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Fig. 1. (a) The phase space for gluon emission in e + e~-ann. is a triangular region 
in the (y, k = lnfc^)-plane. The height of the triangle is given by L — Ins. (b) 
When one gluon is emitted at (yi, Ki) the phase space for a second (softer) gluon 
is represented by the area of this folded surface, (c) Each emitted gluon increases 
the phase space for softer gluons. The total phase space is represented by this 
multifaceted surface. 

r ~ l/k±. Thus the ordering of the gluons is determined by their transverse 
momenta, when measured locally in the emitting dipole rest frame. 

This result can be generalized so that the emission of many gluons can 
be described as a dipole cascade [I]. The phase space for the emissions 
can be represented by the (y,K = ln/c^) diagram shown in Fig. [TJ This 
formulation of the timelike parton cascade is implemented in the Ariadne 
MC [5], which very successfully reproduces experimental data from LEP 
and other e + e~ colliders. 

The phase space in Fig.Q]apparently has a fractal structure. It is possible 
to define a fractal dimension given by D = yj2N c a s /'K [B]. As a s (k"^_) is 
running this is a so called multifractal. It has been discussed if this feature 
is responsible for the "intermittency" signal observed in experimental data. 
However, later it has been realized that a large fraction of the observed 
effect is related to BE correlations. 



Quark confinement can be understood if the colour field is compressed 
to a flux tube by a gluon condensate in vacuum. In an e + e _ -ann. event a 
stringlike field is stretched out between a quark and an antiquark, and when 
enough energy is stored in the field, it can break due by the production of 
a new qq pair. A space-time picture of this process is shown in Fig. [2 

In the Lund model the probability for a definite final state is given by 
the product of a phase space factor and the exponent of a constant times 



3. Hadronization 
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Fig. 2. The hadronization of a high energy (go, <?o) system in a (2, t) diagram. The 
hadrons can be ordered in "rank", 1, 2, ... j, ... n. This ordering agrees on 
average, but not in every case, with the ordering in rapidity. 



the space-time area spanned by the string before it breaks, denoted by A in 
Fig. [2] [7J. This expression can be interpreted as a Wilson loop integral, or 
an imaginary part of the string action. An important feature of the result 
is boost invariance, which is also a property of a homogenous longitudinal 
electric field. 

A gluon carries colour and anticolour charges, and in the Lund string 
hadronization model it behaves as a transverse excitation on the stringlike 
field, stretched between a quark and an antiquark [8j. In a three-jet event 
the string gets a transverse boost, and in the break-up the hadrons are 
produced around two hyperbolae in momentum space, as shown in Fig. 
[9]. Thus fewer particles are produced in the angular region opposite to the 
gluon jet, and this asymmetry was experimentally confirmed, first by the 
JADE detector at the Petra collider [TO]. 

Gluon radiation is singular for soft and collinear emissions. A very 
important feature of the string hadronization model is that it is infrared 
stable. The motion of a soft transverse gluon is soon stopped by the tension 
in the attached strings. In the subsequent string motion the gluon kink is 
split into two corners, which do not carry energy or momentum and which 
are connected by a straight string piece, as shown in Fig. ®l. The energy 
in the small sections close to the quark and the antiquark is not sufficient 
for a hadron, and all breakups will occur in the central string piece, which 
is stretched and breaks up in the same way as the straight string in Fig. [2j 
The string motion with a collinear gluon is shown in Fig. and also here 
the effects of the gluon goes to zero in the collinear limit. 

The situation in Fig. [3] can be directly generalized to many gluons. 
The string is here stretched from the quark to the antiquark via the colour- 
ordered gluons, as shown in Fig. Hb. 
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Fig. 3. The space-time development of a quark-antiquark-gluon event. The string is 
stretched from the quark to the antiquark via the gluon, which moves like a pointlike 
kink carrying energy and momentum. The string breaks by the production of new 
qq pairs, and the final state contains three jets. Soft particles formed in between 
the jets get a boost by the transverse motion of the string. 




Fig. 4. (a): A soft transverse gluon will soon lose its energy. The kink on the string 
is split in two corners and a straight string piece is stretched in a way similar to 
a one-dimensional string, (b): Also for a collinear gluon the energy in the string 
between the quark and the gluon is too small for a breakup of the string, (c): In a 
state with many gluons the string is stretched from the quark to the antiquark via 
the colour-ordered gluons, in the figure from red to antired, from blue to antiblue 
etc. 



4. Spacelike cascades 

As discussed in the introduction, DIS and hadronic collisions are more 
complicated than e + e~-ann.. There are two separate scales, Q 2 and s, and 
two separate problems: inclusive cross sections and final state properties. 
The ladder leading up to the hard interaction (solid lines in Fig. [5]) represents 
the increased parton density in the initial state, and determines the inclusive 
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Fig. 5. A DIS event with ISR, solid lines, and FSR, dotted lines. Virtual links are 
denoted fc, and real emissions qi. 

total and elastic cross sections. The parton links, ki, in these cascades have 
spacelike momenta, and only those branches in the cascade, which interact 
with the projectile, can come on shell and produce real final state particles. 
For exclusive final states also final state radiation has to be added (dashed 
lines in Fig. [5]). This phase is more similar to the cascades in e + e~-ann., 
with timelike virtualities and conservation of probability. 

For high Q 2 and not too small x, the ladder is described by ordered 
DGLAP evolution, where k±i > /cj_j_i and the vertices are determined by 
the quark and gluon splitting functions. For small x, gluon ladders are 
most important, and the evolution dominated by the 1/z pole in the gluon 
splitting function. This pole represents soft emissions, where each step 
in the ladder corresponds to a large step in rapidity. In this section I will 
discuss small x evolution in a semiclassical framework, based on Weizsacker- 
Williams method of virtual quanta. At high energies more than one parton 
in the projectile or the target may interact. This problem will be discussed 
in Sec. [5j 

4-1- Weizsacker-Williams method of virtual quanta 

A Coulomb field, which is boosted to high velocity, is contracted to a 
flat pancake with a dominantly transverse electric field. The pulse will be 
very short in time, and can be approximated by a ^-function: 



E ± ~^ff(t), 



(14) 
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Fig. 6. A Coulomb field in a boosted frame is compressed to a flat pancake 

Here r is the (two-dimensional) distance between the position of the central 
charge and the point of observation (see Fig. E]). The frequency distribution 
is given by the Fourier transform, and consequently approximately constant 
as a function of u, E±(u) ~ g/r. 

The electric field is also associated with an orthogonal transverse mag- 
netic field, with the same magnitude. The energy density in the pulse is 
therefore given by 



The density of photons, or gluons, seen by an observer at point r, is obtained 
by dividing by the energy of a photon, and thus given by 



In the last expression we used that the (twodimensional) Fourier transform 
of a wavefunction proportional to 1/r is given by l/q±- 



A proton is colour neutral, and the colour field from a parton is always 
screened by a corresponding anticharge. Let us study the field from a colour 
dipole formed by a charge at x and an anticharge at y in the transverse 
plane. The transverse field from these charges in point z is given by (cf 




(15) 




(16) 



4-2. Dipoles in spacelike cascades 



Eq. CEH) and Fig. [7Ja) 




(17) 




(a) (b) 



Fig. 7. (a) The transverse colour-electric field in a colour dipole. (b) Gluon emis- 
sion splits the dipole into two dipoles. Repeated emissions give a cascade, which 
produces a chain of dipoles. 



where ri = z — x and r2 = z — y. Denning Y = In a; and R = ri — r 2 , we 
find in analogy with Eq. (|16p the gluon density in point z: 

dn o /ri r 2 \ 2 R 2 , , 



We note that for small t\ we have -R 2 ~ r|, and Eq. (|18j) corresponds to a 
pure Coulomb field oc 1/rf from the charge in x, while at larger distances, 
where R <C T\ « r 2 , the field is screened and falls off ~ 1/r 4 . 

The essential difference between QCD and QED is that the emitted 
gluon carries away colour charge. Thus, if e.g. the gluon with colour rb is 
emitted from an originally rf dipole, the originally red charge is changed 
to blue, and the dipole is changed to a system of two dipoles, a bb dipole 
formed by the originally red (but now blue) charge and the antiblue charge 
in the emitted gluon, and a rf dipole between this gluon and the original 
anticharge. In the large N c limit these dipoles can emit softer gluons in- 
dependently. The number of dipoles increase as a cascade to smaller and 
smaller rapidities Y, as indicated in Fig. [7)3. 

We note that the density proportional to dui/oj = dY corresponds to the 
1/z pole in the q — > qg and g — )• gg splitting functions, which dominate the 
parton distribution for very small x. 



4-3. Double Leading Log approximation 

As mentioned the gluon emission is suppressed when r\ and r 2 are larger 
than R, and the distribution in Eq. (|18p can be separated in a way very sim- 
ilar to the angular ordering in the timelike cascade. We split the expression 
in Eq. (|18j) in the same way as in Eq. ()12|) : 

= l[ Xl + X 2 }. (19) 



R 2 



R 2 



rf + ri 



rr r 



+ (!«-► 2) 



1 '2 



Here the first term in the parenthesis (Xi) is non-singular when r 2 — > 
(and r\ — > R 2 ). Averaging this term over the azimuth angle, (ft, around x, 
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keeping n fixed, we get 

-!- (x 1 d(t>= \e{R- ri ). (20) 

2tt J r{ 
Thus the gluon emission in Eq. (|18|) can be approximated by 
(in a 



d 2 r dY 2vr 



^(i?-r l) + ^?^-r 2 ) 

"l '2 



(21) 



where we have included the proper numerical factor, and used the notation 
a = N c a s /ir. This corresponds to the independent emission from two single 
charges, confined within the regions ri < R. Thus the dipoles are ordered 
in size; the daughter dipole is smaller than her parent. 

A probe with resolution Q 2 can "see" dipoles in a target with size 
r > 1/Q, while smaller dipoles are not resolved. Although non-ordered 
dipoles are not totally excluded in the exact expression in Eq. (fTKj) . the 
approximation in Eq. (j2~T|) implies that the dipoles in a typical cascade be- 
come smaller and smaller. Therefore, for large Q 2 ordered emission chains 
dominate, in which 1/Q < . . . < r j < fj-i < . . . < R, where R is the size of 
the initial dipole in the cascade. Calculating the density of dipoles with size 
r at rapidity y in such a cascade, we get first a contribution from emissions 
directly from the original dipole: 

d v 2 

direct contribution : a — ^- dy (22) 

A two-step contribution is obtained by first emitting a dipole with size r\ at 
rapidity yx, which then emits the observed dipole at a lower rapidity y < y\: 

df 2 pR dr 2 dr 2 
2 steps : a — 7^-dy / a — J- / dy\ = a — ^-dy[a \n(R?lr 2 ){Y — y)\ (23) 
r z Jr rf J y r 1 

Calling the square parenthesis X, we get in three steps 

■ dr2 f R ^ dr i f ¥ f 1 ^ dr 2 f yi j. „ dr2 j, 1 V 2 



3 steps : a^-dy / a— / dyi I a—^- / dy 2 = a^-dy • —X . 

T Jr f\ Jy J r T*2 " V ^ 

(24) 

Summing contributions from n steps, with n = 1 . . . oo, gives the initial 
distribution from a single step with the extra factor 

E^ n = ^o(2v / A). (25) 
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Here Jo is a Bessel function, which for large arguments grows like an expo- 
nential. With Y — y = ln(l/x) we therefore get the result 

~ exp(2 v / aln( J RVr 2 )ln(l/x)). (26) 

This result represents the Double Leading Log, or DLL, approximation, valid 
for small r (meaning large Q 2 ) and small x, when both logarithms are large. 
(We have here assumed a constant coupling a. For a running coupling 
oc l/ln(l/AV), ln(i? 2 /r 2 ) is replaced by ln[ln(l/A 2 r 2 )/ ln(l/A 2 i? 2 )].) 



4-4- BFKL evolution 

When x is small but Q 2 not large, the integral over the ordered dipoles, 
rj, which leads to the factor (ln(i? 2 /r 2 )) n /n! in Eq. (|25|) . becomes small for 
large n. Although suppressed, unordered dipole chains become important, 
in which some rj may be larger than r^—i. We must then use the full 
expression for dipole emission in Eqs. (|18U19p . instead of the approximation 
in Eq. (1211) ■ 

If the rapidity interval Y — y = ln(l/x) is increased by an amount 8Y, 
the density of dipoles J-(Y, r 2 ) will change in the following way: J- can 
increase if a dipole with size r' splits forming a dipole r within the interval 
5Y (a gain term), and it can decrease if a dipole of size r splits into two 
new dipoles (a loss term). This gives the following differential equation: 



dF(Y,r 2 ) a ( [ d 2 r'-r' 2 , 2 f d 2 r' 



2 



2\ 



8Y 2vr \J r 2 (r-r') 2:F(Y,r ^ ' 2 j ^(r-r') 2 ^'" ; j ' 

(27) 

(The gain term has a factor 2, because when a dipole splits one or the other 
daughter can have size r.) 

This equation is equivalent to the LL BFKL equation, conventionally 
formulated in transverse momentum space [11]. An important feature is 
here that the singularity in the gain term at r' = r is compensated by the 
singularities in the loss term at r' = and r' = r. To see this more clearly, 
we note that the integrand in the loss term is symmetric under the exchange 
r' — > r — r'. We can therefore make the following replacement: 



l l 



(r — r' 



1 



r 



12 + (r - r') 2 (r - r') 2 \r' 2 + (r - r') 2 ] ' 



Inserting this expression in Eq. (|27p we see that the singularity in the gain 
and loss terms exactly cancel when r = r'. Making also the variable trans- 
formation k = r/r 2 , we arrive at a conventional form for the BFKL equation 
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in momentum space. The cancellation of the singularity at r' — r = then 
corresponds to the soft cancellation when k^ — kj_ = in momentum space. 

To understand the qualitative features of BFKL evolution we approx- 
imate the dipole distribution in the gain term by its asymptotic form for 
small and large r: 

r' 2 f ~ 4f for r < r' 



r v ~~~ " , (29) 



r 2 (r — r') 2 1 ~ ^ for r > r' 

This approximation is non-singular when r' — r — > 0. In Eq. (|27|) the singu- 
larity in this point was canceled by the loss term, and in this approximation 
we therefore now only keep the gain term. The result is the equation 

W.4^, fl+ f*JV, ) (30) 



r 2 



dY 

We now make the ansatz 



T{Y,r 2 )^e AY (r z )-^-\ (31) 
which inserted in Eq. fl30R gives 



XT = a 



1 1 
- + 



7 1-7 



T (32) 



This approximation reproduces the qualitative features of the LL BFKL 
equation, with singularities at 7 = and 7 = 1. The right hand side 
has an extreme point for 7 = 0.5, which corresponds to A = 4a. The 
approximation somewhat overestimates the contribution from the region 
r' ~ r, and therefore the A- value is larger than the true value A = 4 In 2 a. 
The solution corresponds to an exponential growth for large Y, ~ e XY ~ 
l/x x , which is thus faster than the DLL result given by the exponential of 
the square root of Y. 

The BFKL equation describes the density of partons in a cascade, which 
is relevant for inclusive cross sections. Exclusive final states can be calcu- 
lated in the CCFM model |12[ I13j . which reproduces BFKL evolution in 
terms of weights for final states, in which all real gluons are ordered in 
angle and rapidity. This will be further discussed in Sec. 16.21 



5. Multiple interactions and saturation 

5.1. Experimental evidence 

The strong increase in the parton density at high energy implies that 
a single event often contains multiple parton-parton subcollisions. Such 
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Fig. 8. Distribution in azimuth angle between pairs in events with 7/71"°+ 3 jets 
from the CDF coll. [IB] . The shaded (yellow) region shows expectation from double 
parton scattering. 

events have been observed experimentally |14^ [15] [161 117] . As an example 
Fig. 15. II shows results from CDF for events with 3 jets + 7, which can only 
be described including multiple hard subcollisions. Cf also the talk by Rick 
Field at this school [18] 

5.2. Eikonal formalism 

As mentioned above, rescattering and multiple interactions are most 
easily treated in impact parameter space. The result of repeated scattering 
with momenta kj_j is given by a convolution in k^-space, which corresponds 
to a multiplication in b-space. Thus in impact parameter space the multiple 
interactions are described by a product of the S'-matrix elements for the 
individual interactions: 

S{b) = 5i(6)5 2 (6)5 3 (6). (33) 

If the interaction is driven by absorption into inelastic states i, with 
weights 2/j, the optical theorem gives an elastic amplitude given by 



T=l-e~ F , with F = J2fi- 



(34) 
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For a structureless projectile we then find 
da tot /d 2 b = (2T>, 

a el /d 2 b = (T) 2 , (35) 

V\ne\/d 2 b = (1 - e~^ 2/l ) = Otot - CTgi. 

5. 5. Diffractive excitation, Good- Walker formalism 

If the projectile has an internal structure, the mass eigenstates ^ can 
differ from the eigenstates of diffraction $ n , which have eigenvalues T n . 
With the notation = J2n c kn^n (with \&j n = ^i) the elastic amplitude 
is given by (\E'i|T| , I'i) = J2 c i n Tn = {T), while the amplitude for diffractive 
transition to mass eigenstate Vf^ is given by (^fc|T|^i) = J2 n c knT n cin- The 
corresponding cross sections become 

da cl /d 2 b = C£c 2 ln T n ) 2 = (T) 2 (36) 
da diS /d 2 b = 5^(*i|r|* fc )(* fc |T|*i) = (T 2 ). (37) 

k 

The diffractive cross section here includes elastic scattering. Subtracting 
this gives the cross section for diffractive excitation, which is thus deter- 
mined by the fluctuations in the scattering process: 

da diScx = da m - dff A = «T 2 > - (T) 2 )d 2 b. (38) 

5.4- The BK equation and saturation 

Consider scattering of a dipole with charges at transverse coordinates 
x and y against a dense target at rapidity distance Y. The interaction 
probability is called N(x, y, 1"). Study the change in interaction probability 
when Y is changed to Y + 5Y. The probability that the dipole has emitted 
a gluon at point z, within the interval SY, is given by Eq. (|18p . The change 
in interaction probability is therefore given by |19j 

dN{x,y,Y) a f 2 (x - y) 2 

dY 2W (x-z) 2 (z-y) 2 

[iV(x, z, Y) + N(z, y, Y) - iV(x, y, Y) - iV(x, z, Y)N(z, y , Y)} . (39) 

Here the first two terms in the square bracket give the probability for the 
new dipoles to interact, the third term is the reduction because the orig- 
inal dipole has disappeared, and the last term avoids double counting by 
subtracting the probability that both new dipoles interact. This non-linear 
term prevents the interaction probability to grow beyond 1. 
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If we now take the average, and furthermore assume that (N-N) = {N) 2 
(which may be allowed for a sufficiently dense and homogenous target), we 
arrive at the Balitsky-Kovchegov equation |19j . It is obvious that this equa- 
tion has two fixpoints, given by N = and N = 1. The first corresponds 
to the weak interaction limit, where the quadratic term can be neglected. 
The value N = 1 corresponds to the black disk limit, where the interaction 
probability saturates at the unitarity limit. 



6. Dipole cascade models for high energy collisions 



6.1. Mueller's dipole cascade model 

Mueller's model is based on the dipole evolution discussed in Sees. 14.2 
and 14.41 which describes LL BFKL evolution in transverse coordinate space 
[TTJ [201 EI]- When a dipole emits a gluon it splits in two dipoles, which in 
the large N c limit emit softer gluons independently. The result is a gluon 
cascade in form of a dipole chain, as illustrated in Fig. [7b, where the number 
of links grows exponentially with rapidity as discussed in Sec. 14.41 Gluon 
radiation from the colour charge in a parent quark or gluon is screened by the 
accompanying anticharge in the colour dipole, which suppresses emissions 
at large transverse separation. Therefore the dipoles become on average 
smaller and smaller as the cascade proceeds to smaller rapidities. 

When two cascades collide, a pair of dipoles with coordinates (xj, yj) and 
(xj,yj) can interact via gluon exchange with the probability 2/^, where 



or 



log 



(Xj - yj) 2 {yj - xj) ; 
(xj - xj) 2 (yj - yj) 



2\ i2 



(40) 



We note here in particular that the interaction probability goes to zero for 
a small dipole. This implies that the singularity in the production prob- 
ability for small dipoles in Eq. (|18p does not give infinite cross sections. 
We note also that gluon exchange means exchange of colour between the 
two cascades. This implies a reconnection of the dipole chains, as shown 
in Fig. and the formation of dipole chains connecting the projectile and 
target remnants. 

In Mueller's model the constraints from unitarity are satisfied using the 
eikonal formalism. When more than one pair of dipoles interact, colour 
loops are formed, as shown in Fig.llOl This double interaction is an effect of 
saturation, corresponding to the non-linear term in the BK equation (|39j) . 
It is also related to multiple pomeron exchange and pomeron loops in the 
regge formalism. 

In the schematic illustration in Fig. [TUl rapidity is growing along the 
horizontal direction. We here note that if this event was analysed in a 
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Fig. 9. An interaction between a dipole in the projectile and another in the target 
due to gluon exchange gives a recoupling of the dipole chains. 



Fig. 10. Double interaction results in a dipole loop, corresponding to a pomeron 
loop. 

Lorentz frame closer to the target, the dipole loop could lie completely 
within the evolution of the projectile. Thus double interaction in one frame 
can correspond to a colour loop within the evolution, when viewed in a 
different Lorentz frame. Such loops are not included in Mueller's model, 
and are also not taken into account in the BK equation. 



The Lund model [22^ 123"! [2^1 [25] is a generalization of Mueller's model, 
which also includes: 

- NLL BFKL effects 

- Nonlinear effects in the evolution 

- Confinement effects 

It is implemented in a MC called DIPSY, with applications to collisions 
between electrons, protons, and nuclei. An incoming virtual photon is here 
treated as a qq pair, with an initial state wavefunction determined by QED. 
For an incoming proton we make an ansatz in form of an equilateral triangle 
of dipoles, but after evolution the result is rather insensitive to the exact 




6.2. Lund dipole cascade model 
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form of the initial state. 



6.2.1. Beyond LL BFKL effects 

The NLL corrections to BFKL evolution have three major sources [26j : 

Non-singular terms in the splitting function: These terms suppress large 
z-values in the individual parton branchings. Most of this effect is taken 
care of by including energy-momentum conservation. This is effectively 
taken into account by associating a dipole with transverse size r with a 
transverse momentum k± = 1/r, and demanding conservation of the light- 
cone momentum p + in every step in the evolution. This gives an effective 
cutoff for small dipoles. 

Projectile-target symmetry: A parton chain should look the same if gen- 
erated from the target end as from the projectile end. The corresponding 
corrections are also called energy scale terms, and are essentially equivalent 
to the so called consistency constraint }27] . This effect is taken into account 
by conservation of the negative lightcone momentum components, p_. 

The running coupling: Following Ref . [28] , the scale in the running cou- 
pling is taken as the largest transverse momentum in the vertex. 



6.2.2. Nonlinear effects in the evolution 

As mentioned above, multiple interactions produce loops of dipole chains 
corresponding to pomeron loops. Mueller's model includes all loops cut in 
the particular Lorentz frame used for the analysis, but not loops contained 
within the evolution of the individual projectile and target cascades. As for 
dipole scattering the probability for such loops is given by a s , and therefore 
formally colour suppressed compared to dipole splitting, which is propor- 
tional to a = N c a s /ir. These loops are therefore related to the probability 
that two dipoles have the same colour. Two dipoles with the same colour 
form a quadrupole field. Such a field may be better approximated by two 
dipoles formed by the closest colour-anticolour charges. This corresponds 
to a recoupling of the colour dipole chains. The process is illustrated in 
Fig- fTTT and we call it a dipole "swing" . With a weight for the swing which 
favours small dipoles, we obtain an almost frame independent result. The 
number of dipoles in the cascade is not reduced, and the saturation effect is 
a consequence of the smaller interaction probability for the smaller dipoles. 
Thus the number of interacting dipoles is reduced. Counting only these 
"effective" dipoles, the swing can be looked upon as a 2 —)• 1, or in some 
cases 2—7-0, transition. 
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Fig. 11. Two dipoles with the same colour form a colour octet, which may be better 
approximated by dipoles formed by the closet colour-anticolour pairs. This implies 
a recoupling of the dipole chains. 





Fig. 12. Total and elastic cross sections in pp collisions in the dipole cascade model. 



6.2.3. Confinement 



Confinement is also important. A purely perturbative evolution with 
massless gluons violates Froissart's bound |29j . This is avoided by giving 
the gluon an effective mass. 



6.3. Results 
6.3.1. Total and elastic cross sections 



Results for total and elastic pp cross sections are presented in Fig. [T2l 
Corresponding results for total and quasielastic scattering in DIS are shown 
in Fig. [T31 We here see that the the experimental data are very well repro- 
duced by the model. 
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Fig. 13. Left: Total -f*p cross section for combinations of x and Q 2 , presented as 
a function of the scaling parameter r = (Q 2 /Qq)(x/xo) x [50], with Q = lGcV, 
x = 3 • 1CT 4 , and A = 0.29. Right: The cross section for j*p -> 7p (DVCS) for 
W = 82 GeV as function of Q 2 . Data from HI [31]. 



6.3.2. DifFractive excitation 

Diffractive excitation accounts for large fractions of the cross sections in 
DIS and pp collisions. As mentioned in Sec. 15.31 diffractive excitation is in 
the Good- Walker formalism determined by the fluctuations in the scattering 
amplitude, and we note that the BFKL evolution gives large fluctuations in 
the cascade evolution. 

Study the interaction in a frame, where the projectile is evolved a dis- 
tance y\ and the target t/2 = Y — yi, with Y the total rapidity range 
~ lns/(lGeV 2 ). If we here first take the average over the target states, 
we get the amplitude for elastic scattering of the target. Squaring it gives 
the cross section, when the target is scattered elastically. If we after this take 
the average over the projectile states, we obtain the diffractive scattering of 
the projectile, including the elastic scattering. Thus the expression 

{(T)targ)proj 

-<T) 

targ,proj 

(41) 

gives the cross section for single diffractive excitation of the projectile, with 
the excited mass limited to M\ < exp(yi). Varying y\ gives then da jdM\. 
The resulting cross sections for diffractive excitation in DIS and pp collisions 
are shown in Fig. 1144 together with Zeus data for DIS, and an estimate from 
CDF data for pp collisions. 

It is interesting to study the effects of saturation on diffractive excitation 
|32j . Saturation is not very important in DIS, but in pp scattering the Born 
amplitude is large, and therefore the unitarity effects are also large. Fig. [15] 
shows both the Born amplitude and the unitarized amplitude at 2 TeV for 
different impact parameters b. We see that the width of the Born amplitude 
is very large, and without unitarization the fraction of diffractive excitation 
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(b) M x < 8 GeV 




Fig. 14. Left: The ratio of the total diffractive cross section to the total cross 
section in DIS as a function of W , for M x < 8 GeV. Data from Zeus [33] • Right: 

(cut\ 

The single diffractive cross section for Mx < M x in pp collisions at 1.8 TeV. 
The shaded region is an estimate from CDF data [34 . 




would be correspondingly large. (The smooth lines are fits of the form 
AFP e - aF .) 

However, the unitarized amplitude is limited by 1, and the width of the 
distribution, and therefore the diffractive excitation, is very much reduced. 
This result corresponds to the effect of enhanced diagrams in the conven- 
tional triple-regge approach. We note here also that the strong suppression 
from saturation implies that factorization is broken when comparing diffrac- 
tion in pp collisions and DIS [35J . 

The absorption is most important for central collisions, where thus 
diffractive excitation is most strongly suppressed. As shown in Fig. [16j 
the cross section for diffractive excitation is therefore largest in a ring with 
radius b ~ 1 fm ~ 5 GeV -1 , which grows slowly with energy. 
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b 

Fig. 16. Impact parameter distributions for (T) = (d<Jtot/d 2 b)/2, (T) 2 = da e \/d 2 b, 
and Vt = da^is cx/d 2 b in pp collisions at W = 2 TeV. b is in units of GeV -1 . 



6.3.3. Comparison with Multi-Regge Analyses 

It is also interesting to compare the results from the Good- Walker anal- 
ysis with the multi-regge formalism. To this end we study the contribution 
from the bare pomeron, meaning the one-pomeron amplitude without con- 
tributions from saturation, enhanced diagrams or gap survival form factors. 

When s, M^, and s/M^ are all large, pomeron exchange should domi- 
nate. If the pomeron is a simple pole, we expect the following expressions 
for the pp total and diffractive cross sections: 



"ua = /3 2 (0)( S /sor (0) - 1 = S 2 (0)( 5 /»„) e , 

M i^r ^ HmoM) {ik) «• (42) 



Here a(t) = 1 + e + a't is the pomeron trajectory, and (3(t) and gsp(t) are 
the proton-pomeron and triple-pomeron couplings respectively. Comparing 
our result with this expression we find that it indeed reproduces the triple 
pomeron form, with the following parameter values obtained choosing the 
value so = lGeV 2 for the arbitrary scale parameter |32j : 

a(0) = l + e = 1.21, a' = 0.2GeV" 2 , 

/3 2 (o) = i2.6 mb, m=m^p( 1 2 ; 5 1 t 8t y 

g 3P (t) = const. = 0.3 GeV" 1 . (43) 
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6.3.4. Correlations 



We define the double parton distribution, and the impact parameter 
profile F by the relation 

r(*i, x 2 , b; Ql Ql) = D( Xl ,Ql) D(x 2 , Q\) F(b; x u x 2 ,Ql Q 2 2 ), (44) 

where D(x, Q 2 ) is the single parton distribution. This implies that the cross 
section for double hard interactions at midrapidity is given by 



a 



D 



1 



(^)- (1 + 5 AB ) a cS ' 
with the "effective cross section" , a e g , determined by the relation 

a eS =\[ d 2 b(F(b)) 2 ^ 



(45) 



(46) 



F and a c g are often assumed to depend only weakly on Xi and Q 2 . 
The DIPSY MC shows instead that a spike (hotspot) develops for small 
separations b at larger Q 2 , as illustrated in Fig. [T7] [36]- This result implies 
that <r e ff depends strongly on Q 2 for fixed y/s, as illustrated in Table [TJ Part 
of the correlation is due to fluctuations in the cascade, which can also be 
taken into account in the MC. Without fluctuations jFd 2 b should be 1. In 
Table Q] we see that the fluctuations increase JF by about 10%, which thus 
contributes to the correlations given by JF 2 . 

6.3.5. Final states 

In order to generate exclusive final states, obtained when two dipole 
cascades collide, we first have to determine which dipoles interact and be- 
come recoupled in the way shown in Fig. We note here that BFKL 
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Q(, Qi [GeV*], xi, x 2 


cr cff [mb] 


IF 


1.5 TeV, midrapidity 
10 10 0.001 0.001 
10 3 10 3 0.01 0.01 


35.3 
23.1 


1.09 
1.06 


15 TeV, midrapidity 
10 10 0.0001 0.0001 
10 3 10 3 0.001 0.001 
10 5 10 5 0.01 0.01 


40.4 
26.3 
19.6 


1.11 
1.07 
1.03 



Table 1. Summary of results for a c g and corresponding integrals of the correlation 
function F. 

is a stochastic process, and the interactions between different dipole pairs 
are uncorrelated. This implies that the probability for interaction between 
dipoles i and j is given by 1 — e -2 ^ , where fij is determined by Eq. (f4TJj) . 

As mentioned above, the BFKL equation describes the density of partons 
in a cascade, which is relevant for inclusive cross sections. To describe 
exclusive final states it is necessary to take into account colour coherence 
and angular ordering as well as soft radiation. The latter includes also 
contributions from the z = 1 singularity in the gluon splitting function. 
These effects are taken into account in the CCFM formalism [12\ I13J . which 
also reproduces the BFKL result for the inclusive cross section. 

A very schematic picture of a collision between two dipole cascades is 
presented in Fig. [181 Here three dipole pairs interact, forming two dipole 
loops with an additional loop (denoted A) formed within the evolution of the 
left cascade. Non-interacting branches, like B and C, have to be regarded 
as virtual and must be reabsorbed. 

A reformulation of the CCFM model, called the Linked Dipole Chain 
model, was presented in Ref. [37] , Here it was demonstrated that the inclu- 
sive cross section is fully determined by a subset of the gluons in the CCFM 
approach, denoted "A;_i_-changing" gluons. In Fig. fl~9k . we denote the real 
emitted gluons in a ladder q^, and the virtual links k±i. The fc^-changing 
emissions k±_i are either much larger or much smaller than k\_i—\. This also 
means that q±i ~ m&x(k±i, k±i_i). The chain in Fig. [V9~R is shown in the 
triangular phase space diagram in Fig. 119b . The real gluons qi are ordered 
in p + and in p_, and thus also in rapidity or angle. It was also demonstrated 
that to get the full final states, softer emissions have to be added below the 
horizontal lines in Fig. 119b . as final state radiation. This also includes the 
folds sticking out of the plane, which represent the transverse jets formed 
by the gluons q^. 

Thus in order to generate exclusive final states we should go through 
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Fig. 18. Schematic picture of a collision between two dipole cascades. A dipole 
loop within the evolution is denoted A. Non-interacting branches, like B and C 
have to be regarded as virtual and reabsorbed. 




V 96 
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kJ 






1z 

Q3 




(a) 



Ins 
(b) 



Fig. 19. (a): A parton-parton scattering chain. Virtual links arc denoted fc, and 
real emissions qi. In BFKL dynamics the transverse momenta are not ordered, and 
the result should be the same in any Lorentz frame, (b): The same chain in a (y, 
lnq^) plot. Final state radiation is allowed below the horizontal lines, and on the 
folds representing transverse jets. 



the following steps: 

1. Generate cascades for projectile and target 
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Charged particle ij at 900 GeV, track p± > 500 MeV, for N c h > 6 



Charged particle 17 at 7TeV, track p± > 500 MeV, for > 6 




.+♦+++*• 



ATLAS data 

DIPSV 

Pythia8 



Fig. 20. ^-distribution of charged particles at 0.9 and 7 TeV. The solid line shows 
the result from the D1PSY MC, and the dotted line Pythia. Data from the Atlas 
collaboration 138 . 



2. Determine which dipoles interact 

3. Absorb non-interacting chains 

4. Determine final state radiation 

5. Hadronize 

The main problems in this process are due to the large number of small 
dipoles in the cascades. These have low cross sections, and are therefore not 
a big problem for inclusive cross sections. Because small dipoles correspond 
to high transverse momenta, they do, however, have a large effect on the 
properties of the final states. This implies that the result is sensitive to 
details in the treatment of non-interacting dipoles. Our aim is here there- 
fore not to give very precise predictions, but rather to get insight into the 
dynamical features of small x evolution and saturation. 

As a few examples Figs. 1201 and 1211 show comparisons with ATLAS data 
for minimum bias and underlying events at 0.9 and 7 TeV. Fig. 1201 shows the 
^-distribution of charged particles in minimum bias events. The solid line 
shows the result from the DIPSY MC, and the dotted line the result from 
Pythia. We note that the particle density is well reproduced at 0.9 TeV, 
but does not grow fast enough with energy, which is a problem also for other 
MCs which are not tuned individually for each energy. The properties of the 
underlying event is shown in Fig. [2TI which presents the charged multiplicity 
in the "transverse region" , as defined by Rick Field, as a function of the p± 
of a leading charged particle. The model reproduces quite well the increased 
density for higher p± . More comparisons are found in Ref . [25] . 
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Transverse N c j,g density vs. p^ 1 , \/s — 900 GeV Transverse N c hg density vs. p^£ , \/s — 7 TeV 




2 3 4 5 6 7 8 9 10 2 4 6 8 10 12 14 16 18 2( 

p_ (leading track) [GeV] pj_ (leading track) [GeV] 



Fig. 21. N c h in transverse region to p± of leading charged particle. Notation as in 
Fig. EDI Data from the Atlas collaboration [55] , 

6.3.6. Nucleus collisions 

The model can also be applied to reactions involving nuclei, where e.g. 
saturation effects can be studied including a proper geometry. Some early 
results are presented in Ref. [3D] , 

7. Summary 

In these lectures I have discussed interaction cross sections and particle 
production in e + e~-ann., DIS, and high energy hadronic collisions. It in- 
cludes hadronization, initial and final state radiation, small x evolution and 
saturation. 

I have also presented the Lund Dipole Cascade model for high energy 
collisions, which is based on BFKL evolution and saturation. It is an ex- 
tension of Mueller's model, also including 

• important non-leading effects in BFKL 

• saturation within the evolution 

• confinement 

• A MC implementation DIPSY 

The model gives a good description of inclusive pp and ep cross sections 
(including diffraction), and a fair description of exclusive final states (min. 
bias and underlying event). It has fewer tunable parameters than other 
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event generators, and our aim is not to give very precise predictions, but 
rather to get insight into the dynamical features of small x evolution and 
saturation. As examples it is here possible to study effects of correlations, 
fluctuations, and finite transverse size in a way, which is not easy in other 
approaches. 
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